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ABSTRACT 

We present a method which aUows a direct measureraent of the epicycle frequency k in the 
galactic disk, using the large database on open clusters completed by our group. The observed 
velocity vector (amplitude and direction) of the clusters in the galactic plane is derived from 
the catalog data. In the epicycle approximation, this velocity is the sum of the circular velocity, 
described by the galactic rotation curve, and of a residual velocity, which has a direction that 
rotates with the frequency k. If for some reason the clusters are formed with non-random initial 
perturbation velocity direction (measured for instance with respect to the direction of circular 
rotation), then a plot of the orientation angle of the residual velocity as a function of age reveals 
the epicycle frequency. The data analysis confirms that this is the case; due to the non-random 
initial velocities, it is possible to measure k for different galactic radii. Our analysis considers 
that the effect of the arms on the stellar orbits is small (the Galactic potential is mainly axis- 
symmetric); in this sense our results do not depend on any specific model of the spiral structure, 
like the existence of a given number of spiral arms, or on a particular choice of the radius of 
corotation. The values of k provide constraints on the rotation velocity of the disk and on its 
minimum beyond the solar radius; in particular, Vq is found to be 226 ± 15 kms~^ even if the 
short scale (Rq — 7.5 kpc) of the galaxy is adopted. The mesured k at the solar radius is 42±4 
kms^^kpc^^. 

Subject headings: Galaxy: spiral arms: Galaxy - rotation curve; Open Clusters 



1. Introduction 

Open clusters are ideal objects to study the dy- 
namics of the spiral structure of the Galaxy. A 
large number of them have known distances and 
space velocities (proper motions and radial veloc- 
ities), that are more precisely determined than 
those of individual stars. Another important pa- 
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rameter is their age, which is obtained by means 
of the HR diagrams. In a previous work, Dias 
& Lepine (2005, hereafter DL) integrated back- 
wards the orbits of a sample of open clusters to 
find their birthplaces as a function of time, and ob- 
tained the velocity of the spiral arms. In a similar 
manner, by integrating the orbits, it is possible to 
determine the velocity (direction and amplitude) 
of the clusters at the time of their birth. Since 
the gas of the disk is known to present systematic 
perturbations with respect to circular rotation, or 
"streaming motions" related to the spiral struc- 
ture, it would not be surprising if the new born 
clusters also present some systematic deviations 
from circular motion. The study of initial veloc- 
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ities allows us to verify if the dominant process 
of star formation produces clusters with random 
velocities or, on the contrary, with velocities ori- 
ented in some preferential direction(s). The an- 
swer to this question will certainly be able to put 
restrictions to the mechanisms that trigger star 
formation. 

One result of our investigation is the discovery 
that part of the open clusters have initial veloc- 
ities in a constant narrow angle for long period 
of times. This property allowed us to develop a 
method for a direct determination of the epicycle 
frequency. According to the perturbation theory, 
if an open cluster was born with an initial velocity 
different from that of the circular velocity given 
by the rotation curve, it will follow an orbit that 
can be described by the epicycle approximation. 
The motion can be seen as the sum of the circular 
motion around the galactic center, plus a small 
amplitude rotation around the equilibrium point 
in the circular orbit. The epicycle frequency k is 
a function of the galactic rotation velocity f2 and 
of its derivative <Kl/dR only: 

Based on this equation, the epicycle frequency is 
easily derived from any observed rotation curve. 
However, since there are uncertainties concerning 
the exact shape of the rotation curve, as several 
different curves are proposed in the literature, a di- 
rect measurement of the epicycle frequency, with- 
out passing through equation (1), is of great in- 
terest. Such a measurement, in turn, can help to 
select the best galactic parameters and rotation 
curve. We discuss in this work, in addition to the 
statistics of initial velocities of open clusters, the 
restrictions on the galactic parameters that can be 
derived from the measured epicycle frequency. 

2. The Open Clusters Catalog 

We make use of the New Catalogue of Optically 
visible Open Clusters and Candidates published by 
Dias et al. (2002) and updated by Dias et al. 
(2006) B This catalog updates the previous cat- 
alogs of Lynga (1987) and of MermiUiod (1995). 
The present 2.7 version of the catalog contains 
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1756 objects, of which 850 have published ages 
and distances, 889 have published proper motions 
(most of them determined by our group, Dias et al. 
2006 and references therein) and 359 have radial 
velocities. Recently, Pauzen and Nepotil (2006) 
compared statistically our catalog with averaged 
data from the literature and reached a result which 
is in favor of the data from the catalog. 

3. The method 

The main result of the perturbation theory of 
stellar orbits in a central potential, applied to the 
galactic disk, is that the motion of a star (or of an 
open cluster) can be described as the sum of a cir- 
cular rotation around the galactic center, with an- 
gular velocity f2(i?), plus an epicycle perturbation, 
as illustrated in Figure 1 ( eg. Binney & Tremaine, 
1987, hereafter BT87). The epicycle perturbation 
is an oscillation with frequency k around the un- 
perturbed circular orbit, in both the radial direc- 
tion and the direction of circular rotation: 

r{t) = ra+C{t), with i^{t) ^ b sin{Hit + cj)) (2) 

eft) = 6'o+f^<+ — , with r](t) = a cosUt+cj)) (3) 

(^ and rj represent the displacements in the two di- 
rections, b and a the amplitudes in units of length, 
and (j) the phase). In the frame of reference of 
the guiding center of the epicycle motion, rotating 
with angular velocity 17, the perturbed displace- 
ment describes an ellipse, with the ratio of the 
amplitudes a/b = 2fl/K (BT87), which is approx- 
imately 1.5. Like in the case of a circular motion, 
the perturbed velocity rotates with the same fre- 
quency of the perturbed displacement. The veloc- 
ity components are: 

^' = ^ = 6 K cos(Kt + (j)) (4) 
dt 

, drj 

77 = — = —a K sin[Kt + <p) (5) 

The ratio of the velocity amplitudes is the same 
of the displacement amplitudes. The angle 7 be- 
tween the perturbed velocity vector and the direc- 
tion of circular rotation, taken as a reference (see 
Figure 2), is: 

c' 

7(<) = ton"^(-^) = tan~'^[{b/a)tan{Kt + if)] (6) 
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We introduced a new symbol, ip, for the phase, 
because we want to adopt the direction of circular 
rotation as the origin for 7; the relation between 
the two phases is(^ = (/) + 7r/2. In the figures, 
we represent the galactic rotation in the clockwise 
direction, while 7 increases in the counterclockwise 
direction. The situation 7 = 0° occurs when the 
cluster is at its minimum galactic radius, along the 
epicycle trajectory {(j) = —tt/2, in Equation 2). 

In the hypothetic case a = b, 7 would increase 
linearly with time, with a slope k. In the case a/b 
= 1.5, 7 is a smooth increasing function of time, 
with a slope that on the average is equal to k, 
since in a time T = 2'k/k, 7 varies 360°. If all the 
clusters had a same starting angle, in a plot of 7 
versus age of the clusters the slope of the fitted 
line would give directly the value of k. 

The present day velocity vector of an open clus- 
ter with respect to the Sun, projected on the galac- 
tic plane, is obtained from the data given in the 
catalog: it is the sum of the radial velocity pro- 
jected on the plane, and of the velocity in the plane 
of the sky, derived from proper motion and dis- 
tance, also projected on the galactic plane. We 
first convert the velocities to the Local Standard 
of Rest (LSR). Different sets of solar velocity com- 
ponents {u,v,w) with respect to the LSR are rec- 
ommended in the literature, like the Standard So- 
lar Motion, with {u,v,w) = (-10.4, 14.8, 7.3) (in 
kms-\ Mihalas & Binney, 1981), (-10.4, 5.3, 7.3) 
from the analysis of the Hipparcos data by Dehnen 
et al. (1987), (-9.7, 5.2, 6.7) also based on Hip- 
parcos data by Bienayme (1999), or (-9.9, 12.1, 
7.5) according to Abad et al.(2003), among oth- 
ers. Since differences of a few kms^^ between the 
diverse systems only occur in the v component, we 
varied it to verify how it affects our results. We 
finally adopted v = 9.0 kms~^, since this value cor- 
responds to a minimum of the rms deviation of the 
adjusted points in the determination of the epicy- 
cle frequency. However, this choice is not critical 
for our results; almost the same values of epicycle 
frequencies were obtained with v in the range 6 
to 12 kms^^. Interestingly, we have another argu- 
ment based on the new ideas of the present paper 
in favor of this intermediate value of v between 
the extreme ones, 5 and 15 kms~^. The perturba- 
tion velocities are rotating, and since the clusters 
span a range of ages, even if their initial veloc- 
ities had some preferential direction, the present 
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Fig. 1. — Epicycle orbit around an equilibrium 
point rotating with circular velocity V around the 
Galactic Center. Vp represents the perturbed ve- 
locity. 
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Fig. 2. — The geometry of the velocity vectors. We 
adopt a usual convention with the Y axis oriented 
along the circular rotation at the solar position, 
and the X axis in the galactic radial direction, ori- 
ented towards external regions. One can see that 
^ = e + tan-^{V^/Vy), with 9 = sin-^{y/R) 



direction must be almost random (or cover almost 
uniformly the range -180° to 180°). However, if 
we adopt an artificially incorrect value for v, like 
kms~^ or 20 kms~^, the present day distribution 
becomes either very concentrated around ± 180° 
or around 0°. This experiment tells us that the 
correct v must be near the center of the range 
kms-i- 20 kms-^ 

A further step in the treatment of the velocity 
vectors is their transformation to the local refer- 
ence frame of each cluster, as illustrated in Figure 
2. The new frame has an axis in the local direc- 
tion of rotation and the other in the local radial 
direction. This step includes the corrections for 
differential rotation, described below. The local 
frame of the clusters at present time is the most 
convenient one to start numerical integration of 
orbits. However, for the other use of the veloc- 
ity vectors, which is the direct determination of 
the epicycle frequency, a small correction is still 
required to place the velocity in the frame of ref- 
erence of the guiding center, as later discussed. 

In this work, we do not take into account the 
velocities perpendicular to the plane; only the ve- 
locity components projected on the galactic plane 
are investigated. 

3.1. Corrections for differential rotation 
and local frame of reference 

After electing a rotation curve, we subtract 
from the observed velocity of a cluster the veloc- 
ity expected from pure circular rotation at its posi- 
tion; what is left is the perturbation velocity. This 
is not different from what DL called correction for 
differential rotation. The corrections can be made 
in different ways; the one that we adopted is to 
compute, for the longitude and distance of each 
cluster, the expected U and V velocity compo- 
nents, assuming that the cluster presents pure cir- 
cular rotation, and subtracting the velocity com- 
ponents (0, Vq) of the LSR. Of course, the pertur- 
bation velocity vectors that we derive (observed - 
expected velocities) depend slightly on the choice 
of the galactic parameters Rq, Vq, {dV/dR)Q (the 
radius of the solar orbit, the rotation velocity of 
the LSR around the galactic center, and the slope 
of the rotation curve, respectively). However, we 
are dealing with second order (or differential) ef- 
fects; if we adopt a rotation curve that is too low or 
too high, the LSR velocity and the expected clus- 
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ter velocity arc affected in a similar way, and the 
computed perturbation velocity does not change 
too much. As already discussed, once a perturba- 
tion velocity vector is established, its angle 7 will 
be measured with respect to the local rotation di- 
rection. 

Since the epicycle frequency depends strongly 
on the galactic radius (equation 1), we must select 
the clusters in a narrow range of radius, to be able 
to observe a well defined rotation of the perturba- 
tion velocity. We must make a compromise, since 
too narrow galactic radius ranges will contain only 
a few clusters. After a number of experiments, we 
found that a radial extent of the order of 0.6 - 1 
kpc is convenient. A special radius range is the 
one that contains the Sun (let us say, the range 
i?o± 0.4 kpc), because it is not affected by dif- 
ferential rotation; whatever the rotation curve we 
choose, the expected circular rotation velocity of 
the clusters is the same of LSR rotation around 
the galactic center. 

3.1.1. Rotation curves 

For radius ranges (or rings) relatively distant 
from Rq, the choice of the rotation curve is rel- 
evant. Our task is risky, since we arc willing to 
observe perturbations of the order of 10 kms~^ 
over a line-of-sight projection of rotation velocity 
of the order of 130 kms^^. For instance, if the true 
rotation velocity of the Galaxy in a ring is larger 
than the one that we use for the corrections, the 
perturbation velocities arc imder-correctcd in the 
direction of r]' and the observed 7s will tend to be 
closer to zero than they should be. The effect is 
similar to the one that we described previously, 
concerning the v component of the solar veloc- 
ity. In principle, the fact that under-corrections 
or over-corrections tend to destroy the "normal" 
distribution of present day angles of the velocity 
vectors could be used as a guide to select the cor- 
rect rotation curve. But this is not the purpose of 
the present work. 

Anyway, in order to check the effect of the dif- 
ferential rotation corrections on the derived epicy- 
cle frequency, we performed the calculations over 
a grid of rotation curves. We varied i?o from 7.0 
to 8.5 kpc, in steps of 0.5 kpc, and Vb from 170 to 
250 kms~^, in steps of 10 kms~^. The method is 
very similar to that of DL for the determination 
of the pattern rotation speed; for each point of the 



grid we make use of a rotation curve that passes at 
the corresponding (Rq, Vq) and fits the observed 
CO data taken from Clemens (1985); the same set 
of rotation curves of DL was used here. These are 
curves directly derived from observations, with the 
rotation velocities recomputed for each adopted 
(i?o, Vo). A simple and smooth analytical expres- 
sion was fitted to each of these curves, in order to 
be able to compute the corresponding "theoreti- 
cal" epicycle frequency. We refer to these curves 
as CO-based curves or smooth curves. We also 
performed a number of experiments with constant 
velocity curves, or pure "flat" curves. 

Not all the points of this grid of rotation curves 
and radial ranges produce a valid measurement of 
the epicycle frequency. We later discuss in more 
detail an example with a specific choice of Rq, in 
order to clarify what we consider a valid measure- 
ment. 

3.1.2. The reference frame of the guiding center 

The equations 2 to 5 are valid in the reference 
frame of the guiding center of the epicycle motion, 
which does not coincide with the local reference 
frame at present epoch. A small correction to the 
angle 7, corresponding to a change of the veloc- 
ity component in the rj direction, is required for 
this last frame conversion. The velocity shift only 
affects the direction of circular rotation, since the 
local frame and the the frame of the guiding center 
have the same (zero) radial velocity. 

The expression for the correction is derived in 
Appendix A, based on the equations of motion 
derived by Makarov et al.(2004, hereafter MOT). 
It should be noted that the MOT's equations are 
mainly intended to describe the positions of stars 
as a function of time. In the analysis of positions, a 
secular drift appears and is quite important, this 
being the reason for the stretching of open clus- 
ters. In the present work, we are not dealing with 
positions, but only with velocities (our analysis 
is based on plots of velocity directions as a func- 
tion of age). Velocities are not affected by sec- 
ular terms like positions are; this is why young 
associations that are considerably stretched (like 
Sco-Cen, and the stellar streams) are easily recog- 
nized through their common spatial velocities. In 
our particular case the present day position of the 
particles (clusters) coincide with the origin of the 
individual local frame. As we move towards the 
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past or the future, the local frame rotates with a 
constant circular velocity and the particles move 
out from the origin. 

It is shown in Appendix A that the amplitude 
of the correction 6j to the present day angles, that 
we apply to each point on the graphs like figure 4, 
has an rms value about 15°, and is on the average 
equal to zero. Since the tolerance adopted in our 
analysis (section 3.3) for a point to be considered 
as belonging to a fitted line is 20°, and the fitted 
parameters (k and ip) depend on many points, the 
errors that would be introduced by not perform- 
ing the Sj correction would be small. For such a 
small correction, a first order estimation of it is 
largely sufficient. For this reason, the fact that 
the correction is based on a theory that makes use 
of the Oort's constants, of which we do not know 
a priory the exact values, is not relevant. The use 
of different sets of Oort's constants produces mi- 
nor differences in the corrections (as illustrated in 
Appendix A), that do not affect the fitted param- 
eters. 

3.1.3. The modulus of the velocity perturbation 

The histogram of the modulus of the velocity 
vectors, or residual velocities after correction for 
differential rotation, is presented in Figure 3. The 
peak of the distribution is at about 10 kms~^. In 
parts of the following study we remove from the 
sample the clusters with modulus larger than 20 
kms~^. The reason is that if the amplitude of the 
perturbation is large, the simple epicycle approxi- 
mation is not valid. However, when we are dealing 
with numerical integration of the orbits, which do 
not rely on the epicycle approximation, this re- 
striction is not needed. We also remove from the 
sample the too small velocities (| V |< 3 kms~^ 
), since the velocity components would be of the 
same order of the measurement errors, and the er- 
ror on the angle 7 could be large. 

3.2. The case it!o=7.5 kpc, and clusters 
near galactic radius Rq 

The measurement of the epicycle frequency 
gives us a link between i?o and Vb, like the well 
known link between the same parameters given 
by the Oort's constants A and B. If we fix one of 
these parameters, wc can derive the other. Re- 
cent work often adopts Rq = 7.5 kpc (Racine & 



Harris, 1989, Reid, 1993, and many others). This 
shorter galactic scale, compared to the lAU rec- 
ommended 8.5 kpc scale, is supported by VLBI 
observations of H2O masers associated with the 
Galactic center. The distance to the Galactic cen- 
ter derived from infrared photometry of bulge red 
clump stars (Nishiyama ct al. 2006) is Ro=7.52 
±0.10 kpc, while astrometric and spectroscopic 
observations of the star S2 orbiting the massive 
black hole in the Galactic center taken at the ESO 
VLT (Eisenhauer et al., 2005) gives 7.94 ± 0.42 
kpc. 

Figure 4 (a,b,c) are plots of the angle of the 
perturbed velocity vector with the direction of cir- 
cular rotation, as a function of age. In these ex- 
amples, the galactic positions of the clusters were 
computed considering that Rq = 7.5 kpc, and the 
plots a, b, c, show clusters selected in different 
ranges of galactic radius. In particular, the range 
shown in Figure 4b is 7.1 < R < 7.9 kpc. As we 
already mentioned, this is a special radius range, 
since the choice of the rotation curve has no influ- 
ence on the fitted k. The positions and velocities 
in the galactic plane of the same sample of clus- 
ters are shown in Figure 5. In Figure 4a and 4c, 
the two lines represent the function (6), for a same 
value of K, but different values of ip. When a line 
reaches the maximum angle that can be measured 
(180°), it starts again at -180°. In Figure 4b only 
one line was plotted, as it seems to be a dominant 
one, fitting alone 12 clusters, considering that a 
cluster is fitted if it is situated within ±20", in the 
vertical direction, from the line. The points situ- 
ated around this line are reproduced in Figure 6 
without the folding at multiples of 180°; 2 points 
slightly off the ±20° limit were included. The rms 
deviation of the points from the line is 14°. It can 
be noted from Figure 4 that for ages larger than 
50 Myr, no points considered as not belonging to 
the line arc closer than about 130° (or 9o") from it. 
This justifies our method of "isolating" a line to 
fit the theoretical function. In the present case the 
fitted value of k is 42.0 ± 4 kms^^kpc^^: this is 
the epicycle frequency at the solar radius, one of 
the main results reported here. The initial angle is 
61°, in the frame of the guiding center, or a little 
more, about 80° in the local frame of the clusters. 
This initial angle is in good agreement with what 
is obtained from a totally different method (Figure 
4e) . The probability of the apparent alignment of 
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points being a chance coincidence is discussed in 
the next section. 

The fact that a reasonable number of points 
fall on a same line (for instance, 12 clusters within 
± 20° of the line that begins at 61°) reveals that 
some process confines the direction of the initial 
velocity of the clusters to a same angle with re- 
spect to the direction of circular rotation, for pe- 
riods of time of more than 100 Myr. We postpone 
for the moment the discussion on the nature of 
such process. 

3.3. Probabilities of alignment in a ran- 
dom process 

In the case illustrated in Figures 4b, the sample 
contains 52 clusters. The probability of 12 clus- 
ters being fitted by a theoretical line, as it hap- 
pens in this case, is small. The probability of a 
point falling on a given line which has a "thick- 
ness" of 40° is 40/360 = 1/9. The probabifity of 
12 points being on a given line, out of 52 experi- 
ments, or 12 successes and 40 failures, is [52!/(12! 
40!)](l/9)i2(8/9)4o, « 7x10-3. 

In the example of Figure 4a, with a sample of 
40 points, 2 lines have been fitted, with a total 
of 16 points on them. In a random process, the 
probability of falling on one of the two lines is 2/9. 
If 40 points are randomly placed on the diagram, 
the probability of 16 points falling on any one of 
2 given lines (15 successes and 21 failures) is [40!/ 
(16! 24!)](2/9)i6(7/9)24 « 5 x lO'^. 

In the first example, one could argue that the 

problem cannot be reduced to that of throwing 52 
points on graph and computing the probability of 
12 of them being within a given distance of a line, 
since the line was not previously present; we fit- 
ted their parameters so that they pass through a 
group of points that resulted to be aligned. This 
is equivalent to reduce the number of " successes" , 
for instance, we can choose 2 points to define the 
phase and slope of the line. In this case we require 
only 12 successes and 40 failures; the probability 
is still only about 3%. In reality, the probability 
of the observed alignments of points being spu- 
rious can be discarded, since the fitted value of 
K are always very close to the ones expected from 
theory, and the initial phases are confirmed by an- 
other method. Furthermore, there is a continuous 
variation of k from one ring to the other; for in- 
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Fig. 3. — Histogram of the modulus of the residual 
velocity vectors of the open clusters in the galactic 
plane, after subtracting the normal rotation veloc- 
ity at that radius. 
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Fig. 5. — Observed position and perturbation ve- 
locity vectors of the clusters in the range of galac- 
tic radius 7.1< R < 7.9 kpc. The radius range is 
indicated by full lines and the solar radius by a 
dotted line. In the text wc refer to this region as 
the R = 7.5 kpc ring. The galactic center is not 
contained in the figure. 
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Fig. 4. — Left-side (a,b,c): angle 7 between the perturbed velocity and the direction of circular rotation 
as a function of the age of the clusters. The lines represent the function (6) for different values of ^p. The 
clusters were selected in different ranges of galactic radius, indicated in each frame. The fitted values of k 
were 46, 42 and 45 kms~^kpc~-'- respectively. The right side frames are the histograms of initial velocities 
for the corresponding radius ranges, obtained by numerical integration of the orbits (section 3.4) 



stance, the clusters in figures 4a and 4b are not the 
same, since the galactic rings do not overlap, but 
the fitted slope changes smoothly from one ring to 
the other. We note that the choice of the number 
of fitted lines is arbitrary; similar conclusions can 
be obtained with 1, 2, or 3 lines. We also verified 
that our conclusions do not depend strongly on the 
"thickness" of the lines; if we decide that points 
at distances up to 15° or 25° from a line are fit- 
ted points, the number of fitted points varies, but 
the change in the probability of being on a line by 
chance tends to compensate the variation, in the 
calculation of probability. 

The simple analysis of probabilities tells us that 
we arc not dealing with random alignments, but 
also warns us that we should not attempt to draw 
conclusions from cases with a small number of 
aligned points, like for instance less than half of 
the total number of points being on 3 candidate 
lines, or less than one third of the points being 
on two lines, within a tolerance of 20°. When the 



number of clusters decrease, or when the galactic 
parameters or rotation curves are different from 
the correct ones, we are sometimes confronted 
with situations in which very different values of k 
seem to fit the points with about the same quality 
of the fit. This limits the number of galactic rings 
and the range of parameters that can be explored. 

3.4. Numerical integration of the orbits 
and histograms of initial directions 

The approximation described in the previous 
sections is useful to provide an estimate of the 
epicycle frequency without depending on a precise 
knowledge of the rotation curve. A rough idea of 
the distribution of initial velocity directions is ob- 
tained simply by looking at the starting angles of 
the fitted lines (at age=0). 

To obtain the distribution of initial angles, a 
second method is to perform a numerical integra- 
tion of the stellar orbits. For this integration, a 
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Fig. 6. — Angle 7 between the perturbed velocity 
and the direction of circular rotation as a function 
of the age of the clusters, for the clusters situated 
around the line shown in Figure 4b, with the 180° 
folding removed. The adjusted line gives k = 42.0 
± 4 kins~^kpc~^ and ip = 70°, which corresponds 
to an initial angle 61° (Equation 6). 



first step is to obtain the axis-symmetric potential 
of the Galaxy, by integrating the force V{R)'^/R 
based on the adopted rotation curve. In this way, 
the potential does not depend on any assumption 
on the different components of the Galaxy (disk, 
bulge, etc.). The present day position and veloc- 
ity of a cluster define a quadrivector q(t) contain- 
ing the polar coordinates, the radial velocity and 
the angular momentum, which are the initial con- 
ditions of the integration. In the integration, at 
each time step, q{t + dt) is deduced from q{t) us- 
ing a Runge-Kutta IV method. The value of dt 
is chosen so as to have about 360 points along a 
circle around the Galactic center. The integration 
is performed backward in time, for a total time 
equal to the age. 

If the adopted rotation curve is close to the 
true one, the two methods should produce similar 
initial velocity distributions. Therefore, the com- 
parison between the results of the two methods 
is a tool to make a choice among different rota- 
tion curves. As examples of histogram of initial 
velocity directions, those obtained with the same 
samples used in Figure 4 (a,b,c), are shown in Fig- 
ure 4 (d,e,f, respectively). Note that we expect to 
obtain similar, but not precisely the same distri- 
bution of initial velocities with the two methods, 
since the methods rely on different simplifying as- 
sumptions. For instance, the numerical integra- 
tion takes into account the present day galactic 
radius of individual clusters, while the slope- fitting 
method considers intervals of galactic radius, and 
supposes that the equilibrium radius of the circu- 
lar orbit is approximately the center of the radius 
range. It should be remembered that the initial 
angle in the slope-fitting method is in the frame 
of reference of the guiding centers (Section 3.1.2). 

The agreement between the two methods is 
found to be quite satisfactory when rotation curves 
with high values of Vb (220, 230 kms-^kpc^i) are 
used, but become poor with small Vq, like 180 
kms^^kpc^^. 

3.5. Possible sources of errors 

A source of uncertainty in the value of k, that 
we obtain, in the examples that we discussed, il- 
lustrated in Figure 4, is that k depends on the 
definition of the LSR. It may happen that we are 
not using the best set of solar velocity components 
for the LSR corrections, in the choice described in 
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section 3. The "correct" LSR would be the refer- 
ence frame that rotates around the Galactic center 
with the velocity of the rotation curve at Rq. It 
can be seen from Figure 2 that if we add a cor- 
rection to Vy (in the direction of rotation), the 
computed value of 7 will change. A small change 
in the LSR definition produces a same 5Vy for all 
the clusters, but different angles 7 are affected in a 
different way, so that the slope of the fitted lines in 
Figure 4 may change. The arguments which led us 
to adopt V = 9 kms~^ are explained in section 3.0. 
We verified that typically, if we add 1 kms~^ to the 
Vy component of each cluster, to mimic the effect 
of an error on the v component of the solar mo- 
tion, the fitted 7 is changed about 1 kms~^kpc~^. 
The error on k introduced by the uncertainty in v 
is about ±3 kms~^kpc~^. 

It was mentioned in Section 3 that the ratio 
a/b is approximately 1.5, but the precise value 
depends on the adopted f2 and k. A number of 
tests showed that a small variation of this ratio 
(like for instance taking a/b = 1.4) does not pro- 
duce any change in the fitted values of k, so that it 
was not even necessary to use an iterative proce- 
dure to determine n. Similarly, in section 3.1.2 we 
commented that the correction ^7 for the frame 
of reference of the guiding centers depend on the 
Oort's constants, which in turn depends on the 
adopted rotation curve. In this case too, these are 
second order effects, and an iterative procedure is 
not required. In principle, when a line is fitted to 
the observed points like in Figure 6, the scattering 
of the points due to the errors in radial velocities 
and proper motions contained in the catalog, as 
well as the small effects that we just mentioned, 
are taken into account in the error on the slope 
given by the least-square fitting. We consider that 
our final errors on k are ±4 kms^^kpc^^. 

Other effects like the possible perturbation of 
the orbits by the spiral structure (see eg Lepine 
et al, 2003), or the effects of the local wiggles of 
the rotation curve (section 3.8) are not considered 
here as errors of measurements, but the cause of 
discrepancies between the measured n and theo- 
retical values expected from too smooth rotation 
curves. 



3.6. Detailed analysis and the derived 

value of Vo 

We expect that if a rotation curve is a good 
approximation of the true one, it will produce a 
match between the theoretical and observed k, not 
only near Rq, but at other galactic radii as well. 
In Figure 7 we present k as a function of radius, 
obtained from the same slope-fitting method il- 
lustrated in Figures 4(a,b,c). We remark that all 
the determinations of k rely on this method, the 
numerical integration is used only in the statis- 
tics of initial velocities. The method does not give 
reliable results between 8.0 and 8.9 kpc, because 
the number of clusters becomes too small, and in 
addition, the number of initial phases seem to be 
larger. Only at 8.9 kpc we found again a reliable 
result. Although the observed epicycle frequencies 
do not depart much from those predicted (based 
on Eq. 1) by a smooth rotation curve with Vo= 
220 kms"'^, the slope of function of radius is 

steeper than that of smooth rotation curves. The 
data in Figure 7 was also fitted by a smooth Vq 
= 230 kms"'^ rotation curve modified by the ad- 
dition of a Gaussian minimum centered at 7.8 kpc 
and peak value 6 kms~^ , the corresponding value 
of Vq being about 225 kms"-'^. The existence of 
a minimum is further discussed in a next section. 
We do not intend to present a precise fitting of 
it, but only to illustrate that it could be an ex- 
planation for the anomalous slope of k. Such a 
minimum is almost not perceptible in the rotation 
curve (Figure 10). 

Interesting results are obtained by varying the 
rotation curves with a fixed Rq. Figure 8 presents 
the theoretical and observed epicycle frequency as 
a function of the adopted Vq, for a fixed Ro= 7.5 
kpc. What we call here the theoretical k, is the 
value obtained from direct use of equation (1), 
with f2(i?) given by the adopted rotation curve; 
the observed k is obtained from the slopes of lines 
in plots similar to those shown in Figure 4. For 
each Vo, the same rotation curve is used to deter- 
mine the theoretical frequency and for the differ- 
ential rotation corrections. Figure 8 turns clear 
how the best value of Vq can be derived, and how 
the uncertainty on k affects this Vq value. 

Let us first analyse the case of pure flat rota- 
tion curves (V(i?)=V(i?o)= constant). For such 
curves, it can be derived from equation (1) that 
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K = ^/2Vt. At Rq, VL = Vq/Rq; the dashed hne in 
Figure 8 gives k as a function of Vb, for flat curves 
with Ro = 7.5 kpc. This hne cuts the hne of ob- 
served K at 42 kms~-'^kpc"-'^ and Vq = 226 kms"-'^. 
In other words, if the rotation curve were a flat 
one, our measurement of k = 42 ± 4 kms~^kpc~^ 
would imply that Vq is 226 ±15 kms"""^. 

Besides the pure flat (slope zero) curves, an- 
other family of curves that we consider is that of 
smooth curves, obtained by fitting the CO data 
with different values of Vq (section 3.1.1). The 
line representing the theoretical k for this family 
of curves also cuts the observed line at 226 kms~^. 
The coincidence of the two " theoretical" lines in- 
tersecting at 226 kms"-^ is not surprising, since 
the smooth curves for Vq about 220-230 kms~^ 
are very close to flat ones (see Figure 10). 

3.7. Rotation curves with Rq = 7.5 kpc at 
different galactic radii 

We present in Figure 9 the observed values of 
K (based on figures similar to Figure 4) and the- 
oretical values of k (derived from the family of 
CO-based curves) for two distinct ranges of ra- 
dius, one smaller than Rq, 5.9 < R < 6.5 kpc 
(dashed lines), and the other larger than Rq, 8.2 
< i? < 9.8 kpc (dotted-dashed lines). The lines 
with positive slope are the theoretical ones. The 
two lines corresponding to R ~ 6.2 kpc intersect 
at about Vq = 220 kms~^; while the two lines cor- 
responding to R = 9.0 kpc intersect at about 235 
kms""'^; both values have errors ± 10 kms^^ and 
are in agreement with Vq = 226 kms~^, previously 
obtained from the clusters situated near R = 7.5 
kpc, in Figure 6. 

3.8. Rotation curves with Rq= 7.0 , 8.0 
and 8.5 kpc 

An analysis similar to that of the previous sub- 
sections has been performed, adopting Rq = 7.0 
kpc, Rq= 8.0 kpc and Rq= 8.5 kpc. For each 
adopted Rq, families of smooth rotation curves 
with different values of Vq (as explained in sec- 
tion 3.1.1) were used to derive the theoretical k at 
the radius Rq. On the other hand, the angle of the 
observed velocity perturbation with respect to the 
direction of circular rotation have been plotted, 
using the samples of clusters in the radius range 
Rq ± 0.4 kpc. The study was also made for rings 



Table 1: Values of Vq for which the observed and 

theoretical n are equal, for different adopted Rq 
and different radii of the galactic ring. 



Rq 


ring radius 


best Vq 




kpc 


kpc 


kms~^ 


kms~^kpc— 1 


7.0 


5.7 


195 


50 


7.0 


7.0 


230 


47 


7.5 


6.2 


220 


48 


7.5 


7.5 


225 


42 


7.5 


9.0 


230 


35 


8.0 


6.7 


225 


48 


8.0 


8.0 


230 


40 


8.5 


7.3 


210 


42 


8.5 


8.5 


245 


42 



that do not coincide with Rq] the width of the 
rings was always about ± 0.4 kpc. In all the cases 
the figures look similar to Figure 9, with a rising 
theoretical line and a decreasing line of " observed" 
K. The intersection of the two lines is considered 
as the best estimate of Vq for the given Rq. The 
results are summarized in Table 1. 

It can be seen that the observed n at i?(ring)= 
Rq is almost the same for the different adopted 
Rq (the average of the R= Rq rings in Table 1 
is K=43 ). This is an expected result since the 
samples of clusters are the same in these cases, as 
the positions of the clusters are given relative to 
the Sun. And, as we already mentioned, when we 
study a ring at Rq, the rotation curve has only a 
second order effect in the determination of k. The 
fact that a same best estimate of Vq is obtained for 
Rq= 7.0, 7.5 and 8.0 kpc is an indication that this 
value is a robust result. Furthermore, the fact that 
for Rq = 7.5 and 8.0 kpc almost the same value 
of Vq is obtained from two different rings is an 
indication that the smooth rotation curves that 
we are using in these cases are close to the real 
ones. 

3.9. A wiggled rotation curve and the 
Oort's constcints 

It is well accepted that the rotation curve 
presents irregularities in the solar neighborhood. 
According to Oiling and Merrifield (1998), there 
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is a minimum at about Rq + 3 kpc, and a step 
to higher rotation velocities at about i?o + 4 kpc. 
Honma and Sofue 1997) used the same method of 
Merrificld (1992), and both groups obtain a stcp- 
hke increase of rotation velocity, but at different 
radius: about 1.2 Rq (HS) or about 1.5 Rq (M). 
Brand & Blitz (1993) found the deepest minimum 
at about 1.15 Rq, and the start of a step- like in- 
crease at about 1.3 i?o. Amaral et al. (1996), 
based on an analysis of OH/IR stars, found the 
minimum at 1.07 Rq. We propose the curve pre- 
sented in Figure 10, which contains a minor wiggle 
represented by a Gaussian minimum centered at 
7.8 kpc with a half-width of Ikpc, and depth 8 
kms_i. Such a wiggle produces a reasonable fit of 
the epicycle frequency (sec Figure 7) , and the gen- 
eral curve fits the Clemens (1985) data, excepting 
for the minimum at about 8.5 kpc. It should be 
noted that Clemens did not use the same VLSR 
transformation that we adopted in the present 
work, and that the we applied corrections for the 
different Rq that are only strictly valid for the 
data inside the solar circle. 

In most of the descriptions that we mentioned, 
the minimum is at 1.1 Rq or beyond. So up to Rq 
the curve is still quite flat, and we believe that our 
analysis of the sample of clusters near Rq is not 
seriously affected by the wiggles. Indeed, when 
we plot the angle 7 as a function of age, like in 
Figure 4, but using our curve with aminor wig- 
gle to perform the differential rotation correction, 
the result looks very similar to Figure 4, except 
that the best fit of the slope of the lines gives us 
41 kms~^kpc~^, instead of 42 kms~^kpc~^. This 
difference is well within the error bars. 

An important aspect of the minimum in the 
curve is that it could explain the apparent con- 
tradiction between our results and previous re- 
sults based on the Oort's constants. The value 
of Vo that we obtain may seem too high, since 
it corresponds to Vo/Rq = 30.1 (the units are 
kms~^kpc~^ in what follows). In a first order ap- 
proximation, the parameters Vq and Rq are linked 
together through the Oort's constants A and B, 
which are determined by observations (Vq/Ro = 
A - B, (dV/dR)7^o = -A-B ). A value of A-B which 
seemed to be well established is 26, confirmed by 
recent observations of a different nature (Kalirai 
et al., 2004) which give 25.3. Feast & Whitelock 
(1997) obtained 27.2 based on Hipparcos proper 



motions of Cepheids. However, Oiling and Mer- 
rifield (1998) verified that the Oort's constants A 
and B differ significantly from the general Vq/Ro 
dependence, in the solar neighborhood. They at- 
tributed this effect to an anomaly in the local gas 
distribution. Oiling & Dehnen (2004) showed that 
the most reliable tracers of the "true" Oort's con- 
stants are the red giants, and derived A-B = 33. 
Another previous determination with high value of 
A-B is that of Backer & Sramek (1999), who mea- 
sured the proper motion of Sgr A, at the Galactic 
center, equal to 6.18 ± 0.19 mas yr~^, equivalent 
to 29.2 ± 0.9. This is an interesting result since it 
does not depend on local wiggles. Branham (2002) 
obtains fl= 30.3, Fernandez et al. (2001),Q =30, 
and Miyamoto and Zu (1998), Q = 31.5, all of 
them from the kinematics of OB stars. Meztger et 
al. (1998) obtained 31 from Cepheid kinematics, 
and Mendez et al. (1999) obtained 31.7 from the 
Southern Proper Motion Program. We conclude 
from this discussion that the value of Vq that we 
obtain, and our choice Rq = 7.5 kpc, are not in 
contradiction with the constraints of the Oort's 
constants. 

As we already mentioned, one of the main prop- 
erties of the wiggles is a minimum at about 1.1-1.2 
Rq, observed by different authors and using differ- 
ent tracers. We only insist here on the existence 
of the minimum because it has not always been 
recognized, possibly because there seems to be no 
obvious reason for it. The existence of this sharp 
(half-width about 1.0 kpc) minimum could possi- 
bly be related to corotation. This resonance is a 
frontier between two regions of the galactic disk, 
one in which the matter flows towards the cen- 
ter, and one in which the matter flows towards 
external regions. This tends to form a "vacuum" 
around corotation, as confirmed by hydrodynamic 
simulations (eg. Lepine et al. 2001). If less gas 
is present, less stars form, and after a few billion 
year a "gap" ring with a smaller stellar density 
forms. 

3.10. Statistics of initial directions 

Typically, a ring with a width about 0.8 kpc 
contains 3 or 4 main directions of initial velocity 
perturbation, seen as peaks in the histograms. A 
preliminary analysis of the position on the galac- 
tic plane where open clusters were born, indicates 
that at least two different directions can originate 
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from a same position in the galactic plane, at a 
same epoch. It was shown by DL that the clusters 
were born in spiral arms. We can now say that a 
same spiral arm can prodncc star clusters going in 
different directions. However, this kind of analysis 
is difficult to perform, since the samples become 
very small when narrow bins of ages and initial 
positions are considered. 

Different rings may have different main direc- 
tions, but some general tendencies can be ob- 
served. Figures 11 to 13 present the distribution 
of initial velocities in the local (U,V) frame of each 
cluster, for different ranges of radius, obtained by 
numerical integration. One is the local region, 
from 7.1 to 7.9 kpc (Figure 11), which has already 
boon studied in detail, and is supposed to include 
the corotation radius. Another region, from 5.0 
to 7.1 kpc (Figure 12) is inside the corotation cir- 
cle. Wc consider a large galactic radius interval 
to obtain better statistics. Finally, in Figure 13 
we show the distribution for clusters in the radius 
range 8.1 to 12 kpc, external to the corotation cir- 
cle. One can see some preferential or non-random 
directions, which are the basis of the present work, 
since they offer the possibility of measuring the 
epicycle frequencies. Although the statistics are 
poor, it does not seem that a description in terms 
of pure ellipsoids would be an adequate one. 

Clear differences appear between the three re- 
gions, with dominant positive U velocities at large 
galactic radii, and dominant negative U velocities 
in the internal region. There seems to be a change 
of about 180° in the main velocity direction in the 
comparison of internal and external regions (from 
-50° to + 130°). This is in general agreement with 
the idea that outside corotation, the spiral arms 
tend to pull the galactic material towards the ex- 
terior, while inside corotation the arms pull the 
material towards the center. 

In principle, we might expect that the clusters 

start their life traveling along the spiral arms, so 
that they would not leave the arms before the 
brightest stars reach the end of their life. This 
would contribute to the well defined visible aspect 
of the arms. DL showed that if we select clusters 
with ages smaller than 20 Myr, the spiral struc- 
ture appears clearly, while for ages larger than 30 
Myr, it almost disappears. This time scale is com- 
parable to the time required for a velocity vec- 
tor to change its direction by 90°, which is 1/4 



of the epicycle period. Interestingly, the domi- 
nant initial directions for the internal regions of 
the Galaxy satisfy the condition above. If we con- 
sider the vector sum of the circiilar velocity and 
of the perturbation velocity, for instance a vector 
with modulus about 230 kms~^ in the V direction 
plus a perturbation with modulus 15 kms""'^ in the 
direction -70° or in the direction -140° (main di- 
rections seen in Figure 12), the resulting velocity is 
approximately along the spiral arms, considering 
a typical pitch angle of about 14°. 

4. Conclusions 

The study of samples of open clusters situated 
in small intervals of galactic radius or rings with 
width < 1 kpc, reveals the presence of several se- 
quences of 4 to 12 clusters that have been formed 
with a same initial direction of the perturbation 
velocity, as measured with respect to the direction 
of circular rotation. The expected rotation of the 
perturbation velocity vector from the epicycle the- 
ory can be observed in plots of the angle versus age 
of the clusters. In such plots, the sequences of clus- 
ters along lines that correspond to a same initial 
angle span ranges of ages of about 40 to 150 Myr. 
In a given ring, usually 2, 3 or more such preferen- 
tial directions are observed. Probably, the initial 
angles are dictated by the local hydrodynamics in 
the spiral arms. The observations suggest that 
portions of the shock waves associated with the 
arms are coherent, or keep the same characteris- 
tics, during period of times of more than 100 Myr. 
Therefore, the observations provide an argument 
in favor of star formation triggered directly by the 
spiral shocks, and not by secondary processes like 
the explosion of supernovae. 

If we average the initial direction of the per- 
turbation velocity over a relatively large range of 
galactic radius (a few kpc), the individual coherent 
groups become less prominent, or average out, but 
a broad general distribution of initial directions 
appears, with a peak at about -70°, for galactic 
radii smaller than Rq (or smaller than the corota- 
tion radius). This particular initial direction helps 
the arms to keep their sharply defined aspect, since 
the clusters do not move out from the arms before 
the bright short-lived stars reach the end of their 
life. 

The property of the cluster formation mecha- 
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nism of producing initial perturbation velocities in 
a given direction for groups of clusters spanning a 
range of ages allows us to perform a direct deter- 
mination of the epicycle frequency k in the galactic 
disk. This measurement shows only a minor de- 
pendence on the rotation curve which is used to 
perform the corrections for differential rotation, to 
find the perturbation velocity. For a given rota- 
tion curve, there are two ways of obtaining the 
epicycle frequency at a radius R. One is to make 
direct use of the expression for k as a function of 
f2 and its derivative, and the other uses the plots 
of the angle versus age of the open clusters. We 
compared the results obtained using different ro- 
tation curves, based on different values of Rq and 
of V{) . If we fix for instance i?o and vary Vo , we can 
find a best value of Vb, for which the "observed" 
and the "theoretical" values of k coincide. 

If i?o = 7.5 kpc is adopted, using the sample 
of clusters in the range 7.1 < i? < 7.9 kpc, we 
find Vb= 226±15 kms~^. The same value of Vq 
is obtained when we adopt i?o = 7.0 kpc or i?o = 
8.0 kpc. Therefore, Vq= 226±15 kms~^ is a robust 
result, which does not depend on the precise choice 
of i?0) in the 7-8 kpc range. 

We examined samples situated in rings with dif- 
ferent galactic radii; for instance with mean radius 
between 6 and 8 kpc. The rotation curve i?o = 7.5 
kpc , Vo= 226 kms^^, with a minimum at 7.8 kpc 
is found to be a satisfactory one, since it predicts 
values of k that are consistent with the observed 
ones, not only at i?o, but also at these different 
galactic radii. The ratio V0/-R0 that we recom- 
mend is 30.1 kms^^kpc^^, which is larger than 
that "classical" value obtained from the Oor's con- 
stants (A-B), 25 kms^^kpc^^. However, as we 
discussed in a previous section, there are many 
recent measurements that result in Vo/i?o close to 
30 kms~^kpc~^, including a method that does not 
depend on local irregularities of the rotation curve. 
The method to determine Vq that we proposed is 
a new one, which is something much needed in 
view of the conflicting results that have been pub- 
lished in the last decade, almost always based on 
analyses of Oort's constants. 

The work was supported in part by the Sao 
Paulo State agency FAPESP (fellowship 03/12813- 
4) by the Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico. 



5. Appendix A 

The present paper deals with open clusters that 
are spread over distances distances of several kpc, 
so that it is convenient to use an individual frame 
of reference for each cluster, to describe their mo- 
tion. The local frame is chosen to rotate around 
the galactic center with the velocity of the rotation 
curve at the radius where the cluster is at present 
epoch. In such a reference frame, in which ^0 = 
and ?7o =0, the MOT's equations of motion (see 
section 3.1.2) reduce to: 

,A-B . , , 
+% sin{K.t) 

where A and B are the Oort's constants. Taking 
the derivative: 

^' ='^sin{Kt)+£,'^C08{Kt) 



V = -^sm{Kt) + Vo^- Vo Q cos{Kt) 

The only non- harmonic term is ?7o-^, a term 
that does not increase with time and represents a 
shift in velocity, which corresponds to the change 
from the local frame to the frame of the guid- 
ing center (the frame in which the motion is de- 
scribed by pure harmonic functions in both axes). 
At the instant t=0 the velocity in the local frame 
is (^Q,??o) and in the frame of the guiding center 

The initial angle 7 in the frame of reference of 
the guiding center is: 

-/OGC = atan{ % ) 

vol-'- ~ Bl 

SO that: 

57 = atan{C x tan['yo\) — 70 

where S'y is the correction to 70 needed to con- 
vert it from the local reference frame to that of 
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the guiding center. The constant C obtained us- 
ing the classical values of the constants A =15 
kms~^kpc~^ and B= -11 kms~^kpc~^ (Mihalas & 
Binney, 1981, chapter 8) is 0.42; using the Oort's 
constants adopted by MOT, A= 12 knis~^kpc~^ 
and B= -15 knis~^kpc~^, C= 0.55. since only an 
order of magnitude is needed, we can adopt C= 
0.5. The function ^7 is illustrated in Figure Al; 
its average is zero, its maximum about 24° and 
rms value 15° in the worst case (C= 0.42). 
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Fig. 7. — Epicycle frequency measured at different 
galactic radii. The dashed lines represent the the- 
oretical K corresponding to smooth rotation curves 
with Vo = 230 kms"^ (upper curve) and 220 kms~^ 
(lower curve) both with i?o = 7.5 kpc. The dotted 
line correspond to the theoretical k obtained with 
a smooth rotation curve with Vq = 232 kms~^, to 
which a Gaussian minimum centered at 7.8 kpc 
has been added. 




Fig. 8. — The epicycle frequency as a function of 
Vo for different rotation curves, all based on Rq 
= 7.5 kpc. The full line with positive slope is the 
theoretical prediction for the series of CO-based 
rotation curves with different values of Vq (see sec- 
tion 3.1.1). The dashed line is the theoretical k 
for pure flat curves with different Vq. The full line 
with negative slope shows the observed values of 
K, based on angle versus age fitting, like in Figures 
4(a,b,c). The theoretical and observed k intersect 
at at about 226 kms~^. 
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Fig. 9. Theoretical and observed values of k as a 
function of Vq of the CO-based rotation curves (all 
with Ro = 7.5 kpc), for two distinct galactic radii: 
6.2 kpc (dashed lines) and 9.0 kpc (dashed-dotted 
lines). The positive slope lines are the theoretical 
ones, and the negative slope are the observed ones, 
like in Figure 8. 
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Fig. 10. The rotation curve of the Galaxy for Rq 
= 7.5 kpc. The line is a fit to the CO data, using 
the empirical expression V= 240 cxp(-(3.6/r)^) + 
350 exp(-r/3.5 -0.12/r); the units are km~-^ and 
kpc. A very small wiggle was added, represented 
by a Gaussian minimum -8 exp(-[(r-7.8)/1.0]^), in 
order to better fit the epicycle frequccy data (Fig- 
ure 7). Including this term, Vq — 226 kms^^ pre- 
cisely at Rq, but outside the minimum V =230 
kms~^. 
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Fig. 11. — Distribution of initial velocities of clus- 
ters in the galactic radius range 7.1 to 7.9 kpc 
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Fig. 12. — Distribution of initial velocities of clus- 
ters in the galactic radius range 5.0 to 7.1 kpc 
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Fig. 13. — Distribution of initial velocities of clus- 
ters in the galactic radius range 8.1 to 12.0 kpc 
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Fig. 14. — Correction to the angle 7 to move it 
to the reference frame of the guiding center of the 
epicycle motion, as a funtion of 7. The dashed 
line is with C=0.55, full line with C=0.42 (see 
Appendix). 
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